Current direction induced rectification eff'ect on (integer) quantized Hall plateaus 
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Current polarization induced rectification of the quantized ISall plateaus (QHPs) is studied within 
a Hartree type mean field approximation for asymmetrically depleted samples. We first investigate 
the existence of the current carrying incompressible strips (ISs), by solving the self-consistent equa- 
tions, and their influence on magneto-transport (MT) properties. Next, the widths of the ISs are 
examined in terms of the steepness of the confining potential profile considering gate defined Hall 
bars. The corresponding MT coefficients are calculated using a local Ohm's law for a large fixed 
current and are compared for symmetric and asymmetric depleted samples. We predict that, the 
extend of the QHPs strongly depend on the current polarization, in the out of linear response 
regime, when considering asymmetrically depleted samples. Our results, concerning the extend of 
the QHPs depending on the current polarization are in contrast to the ones of the conventional 
theories of the integer quantized Hall effect (IQHE). We propose certain experimental conditions to 
test our theoretical predictions at high mobility, narrow samples. 



PACS numbers: 73.20.-r, 73.50.Jt, 71.70.Di 

Surprisingly, open questions remain even nowadays in 
the theory of the IQHE almost three decades after its 
discovery [T]. When a two dimensional electron system 
(2DES) is subject to a strong perpendicular magnetic B 
field, the energy spectrum is (Landau) quantized. Due 
to the gapped density of states (DOS), the measured 
longitudinal and Hall resistances, and Ru, present 
anomalies if the electron density (ri(,i) is an integer mul- 
tiple of the quantized magnetic flux density (n^), such 
that Rl = and Hall resistance becomes quantized, i.e. 
Rh = where filling factor v{= Uei/n^) is an inte- 
ger, e is the electron charge and h is Planck's constant. 
In a first order approximation, two main schools have 
emerged in giving an explanation to the IQHE, namely 
the bulk [2] and the edge [3j HI |5] pictures, which are 
thought to be unavoidably in contrast to each other in 
answering the question " where does the current flow 
?" . Moreover, there is an ongoing debate about whether 
transport at the edges occur in the compressible [3l IS] or 
in the incompressible [6l|7l|8j|9] states, which are formed 
as a direct consequence of Landau quantization and their 
widths are determined by the Coulomb interaction. The 
local probe experiments present a strong evidence sug- 
gesting that the current is carried by the incompressible 
strips |10|, lllj. In particular, experiments performed at 
the von Klitzing's group where a scanning force micro- 
scope was used to measure the spatial distribution of the 
Hall potential across the 2DES as a function of the B 
fleld [lOj . The observed dependence of the potential pro- 
flle on v already suggests the dominant role of the e — e 
interactions, leading to finite widths of both compress- 
ible and incompressible strips (ISs), where the current 
is carried by the later. The Hall potential profiles were 
categorized mainly to three types: Type I, the potential 
varies linearly in position 1.6 < v < 2, Type II non- 
linear spatial variation w 2) and. Type HI, where the 



potential strongly varies at the edges, however, is con- 
stant at bulk, 2.05 < v < 2.3. The observations were 
explained within the self-consistent (SC) Thomas- Fermi- 
Poisson theory of screening [12 plus the local Ohm's law 
(LOL) which are the bases of the present work to be 
discussed later. In a subsequent theoretical work [8j the 
effect of finite extent of the wave functions on the incom- 
pressible strips (ISs) was simulated by a spatial averaging 
of the local quantities over quantum mechanical length 
scales such as the Fermi wavelength A or magnetic length 
I (= \fTifeB. The spatial averaging also enabled them to 
relax the strict local approximation considering MT and 
to lift the artifacts arising from the Thomas-Fermi ap- 
proximation (TFA) . The main outcome of this work was 
to show explicitly that, if there exists an IS somewhere 
in the sample (which is translation invariant in the cur- 
rent direction) the system is in the QH regime, i.e. the 
widths of the QHPs strongly depend on the widths of 
the ISs. Moreover, its predictions on the asymmetry of 
the QHPs with respect to the classical Hall resistance 
curve, depending on the mobility and sample width, are 
confirmed experimentally |13j . 

In the present work, we first present our geometry and 
the related electrostatic problem, which in turn deter- 
mines the width of the ISs. Next, we investigate the effect 
of a large current on the local electron density nci(x) us- 
ing the LOL. The DOS D(E) and the MT coefficients are 
obtained from the SC Born approximation [14^ . The cur- 
rent density j(r), nci(x) and Rh are compared in the out 
of linear response regime (LRR) for generic and asym- 
metrically depleted samples. The asymmetric distribu- 
tion of the ISs with respect to the center of the sample 
and its effect on Rh is utilized as an experimental test 
for two different boundary conditions, thereby density 
profiles. At asymmetrically depleted samples, in which 
the potential profile is steeper at one side than the other. 
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FIG. 1: (Color online) (a) The spatial distribution of ^{x) 
and the corresponding jy{x) considering a generic (symmetric, 
i.e. hi = br = 150 nm, assuming Vl = Vr = 0) 3 /^m sample, 
(b) The SC and electrochemical potentials under low (solid 
lines) and high bias (broken lines). All calculations are done 
for h{eB/m)/E% (= Q,/E%)=QM, where E% (=12.75 meV) 
is the Fermi energy at the center, i.e. X = 1.5 /im, at the 
default temperature fcsT/n = 0.025. A homogeneous donor 
density of no = 4 x 10^^ cm~^ is assumed, with the impurity 
parameter V /Q, = 0.03. 



we show that the IS at the steep edge is narrower than 
the one at the smoother edge, which in turn determines 
the current distribution and the extend of QHPs. As a 
result of the applied current a Hall potential develops 
within the sample, whose slope is determined by the cur- 
rent direction and amplitude. Therefore, if the Hall po- 
tential added to the electrostatic potential has the same 
slope sign with the SC potential where the narrow IS re- 
sides, this IS is enlarged and the QHP becomes wider 
otherwise, becomes narrower. Such a rectification effect 
due to current polarization is counter intuitive consider- 
ing the conventional theories, both bulk and the edge, 
since the extend of the QHPs mainly depend on the mo- 
bility at fixed temperature and current amplitude, not 
its direction. Here we follow the path of Ref. ^ in de- 



scribing our 2DES considering the historical Chklovskii 
geometry |15j . i.e. a translation invariance in y— (cur- 
rent) direction, where donor distribution no is assumed 
to be homogeneous residing together with the electron 
layer on z = plane and the 2DES is depleted from the 
edges by applying Vl and Vr to the metallic gates on 
sides. In the screening theory of the IQHE, the Coulomb 
interaction is included to a spin degenerate single particle 
Hamiltonian, within the a Hartree type approximation, 
via adding an effective mean field potential given by 
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nci(t)K{x, t)dt, 
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where k is an average dielectric constant (=12.4 for 
GaAs) at the interface of the 2DES and the kernel K{x, t) 
preserves the boundary conditions V{—d) = Vl and 
V{d) — Vr for the above described model. The total 
potential (energy) is then 



V{x) = Vt,g{x) + VG{x) + VH{x) 



(2) 



where Vbg{x) is the background potential generated by 
the donors and Vg{x) by the gates for a sample width 
of 2d. To calculate the Hartree potential one needs the 
electron density distribution, which is given within TEA 

by 



nci(x) = / dED{E) 
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where ks is the Boltzmann constant, T the tempera- 
ture, pl{x) = ^*^^ — V{x) the electrochemical potential and 
/i*q the chemical potential at equilibrium. We define the 
widths r of the broadened Landau levels (LLs) from the 
mobility dependent short range broadening |8l [Ij . The 
SC scheme is closed by the Eqns. |2] and [3] provided that 
the left and right depletion lengths, 6; and bj., are given. 
The numerical task is now to solve these equations by 
iteration until the electron density distribution remains 
unchanged up to a numerical accuracy of 10~*. In the 
next step the local current density j(r) is calculated as- 
suming a fixed current in y direction / = J_j^jy{x,y)dx 
via Ohm's law 



VA.*(r)/e = E(r) = p(r)j(r), 



(4) 



provided that the resistivity tensor p(r) is known through 
the DOS 13 H] and assuming a stationary state using the 
local electric field E(r) obtained in the previous step. 
The translation invariance is utilized together with the 
equation of continuity V • j(r) = and V x E(r) = to 
obtain 



J. = 0, Ey{x) = 

jy{x) = EI/pl{x), E,{x) = E°yPH{x)/pL{x), 



(5) 



where Pl{x) and Ph{x) are the diagonal and off-diagonal 
entries of the resistivity tensor, respectively, and the con- 



stant electric field in y direction Ey = I. 



J-d 



dx 



3 



1 Xft"") 2 3 




X(nm) 



FIG. 2: (Color online) The local variation of the filling factor 
and the current densities for three selected B values, indicated 
in Fig. |3j considering: (i) generic (black solid line) and (ii) 
asymmetric samples, by setting Vl = and Vr = —1.1 V 
which results in b; = br/2 — 75 nm. Current polarizations 
are depicted by (red) broken lines for (+) and (blue) dash- 
dotted lines for (-). The asymmetric depletion is shown by the 
diagonal shaded region on LHS and horizontal shaded region 
on RHS (c), wich results in Ep — 13.12 meV. Insets show the 
regions of ISs on the right side. 



Then fJ-*{x) (now position dependent) is obtained from 
Eqn.|4]by integration, up to a constant which is fixed by 
TT-ci- In our numerical scheme we start with nci{x) calcu- 
lated without current, then calculate the current distri- 
bution for a given fixed /. Next, we obtain ^i*{x) such 
that riei is kept constant and start the new iteration from 
the newly calculated ne\{x). This procedure is continued 
until convergence is obtained. In this paper we apply 
the above described calculation scheme to an asymmet- 
rically depleted gate defined sample. Following Ref. [S], 
we perform a spatial averaging over A (~ 33 nm) to simu- 
late the effects of the finite extend of the wave functions, 
which also lifts the local strictness of the Ohm's law. We 
show that the large current induces an asymmetry on the 
widths of the ISs due to the tilting of the Landau levels 
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FIG. 3: (Color online) The spatial distribution of the ISs 
(blue regions) as a function of B for generic (a) and asym- 
metric samples with (-I-) (b) and (-) DC polarization (c). The 
corresponding Rh and Rl (d). Vertical lines indicate the B 
field values shown in Fig. [2] 



as a result of self-consistency, i.e. adding the Hall poten- 
tial to the total potential and recalculating the electron 
density. The amplitude of the current / ~ 3 /iA is suf- 
ficiently large being in the out of LRR. Our aim is, first 
to present this current induced asymmetry calculated for 
a generic sample, which is equally and largely depleted 
from both edges. Next using the side gates, we deplete 
the sample asymmetrically such that the potential on the 
right hand side (RHS) is smoother than that of the left 
hand side (LHS). Therefore, the IS on RHS is larger even 
without any current induced effect. In the last step we 
apply a large negative (-) and a positive (-I-) DC current 
to the system and investigate its effect on the density 
distribution and the widths of the QHPs. 

The applied current introduces a Hall voltage, which is 
added to the SC potential, hence tilts the LLs and elec- 
trons are redistributed accordingly. As a consequence of a 
(-I-) bias, if there exists an IS, the spatial extend of the en- 
ergy gapped region on the RHS becomes wider, whereas 
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shrinks on the LHS, resulting in a wider IS on the RHS 
and a narrower IS on the LHS. Such a situation is shown 
in Figjl^. The local electron distributions (or equiva- 
lently the local filling factor ^{x) = 2TrP) calculated at 
low (thick solid lines) and high (thick broken lines) cur- 
rent biases is shown together with the current density 
distribution (thin vertical lines) in the upper panel for 
a generic sample. It is clear that, the IS on the RHS is 
larger than the one on the LHS under the large current 
bias (cf. the inset of Fig. [T| and the current (horizontal 
lines) is well confined to the ISs. Fig.[l|D presents the SC 
potentials (thick lines) together with the position depen- 
dent electrochemical potentials (thin lines). We observe 
that at the large bias the Hall potential tilts the LLs. 
Note that, since the compressible regions can almost per- 
fectly screen the Hall potential, the major effect on the 
fi* {x) is observed at the regions where an IS resides. Now, 
we investigate the effect of current bias induced density 
asymmetry at the asymmetrically depleted samples. In 
Figj2] (black) soHd curves present v^x) of a generic sam- 
ple. For the lowest B value (a), no ISs exist larger than 
A (see the inset) for the symmetric and negatively biased 
asymmetric (dash-dotted blue curves) samples. There- 
fore, the electron and the current densities both remain 
symmetric (note that the system is completely compress- 
ible) , hence the induced Hall potential can be almost per- 
fectly screened. The current distribution exhibits, local 
spikes at the positions of ^{x) ~ 2 since, pl{x) assumes 
very small values in the very close vicinity of I'i^x) — 2, 
although no IS exists. One can clearly see that some 
amount of current is still flowing from the bulk. How- 
ever, for the asymmetric sample under (-1-) bias (broken 
red lines) the IS on the RHS is larger than A, hence, the 
current is confined within this region mainly, meanwhile 
no current flows from the bulk. We observe that for this 
B value, the generic and asymmetric (-) samples are out 
of the QHP. When increasing B slightly (Fig.[2]D), an IS 
is now well developed on the RHS of the generic sam- 
ple where most of the current is confined to, however, 
due to the local minima at the LHS mention before some 
current also flows form this side. The interesting point 
is that now no current flows from the bulk (up to our 
numerical accuracy, i.e, 10""'^^ A) and the system is in 
a QHP both for generic and asymmetric (-I-). Increas- 
ing the B field furthermore, results in formation of an 
IS also at the asymmetric (-) sample, where all three 
samples are in the plateau regime. When the center h' 
becomes very close to two (not shown here) the two ISs 
merge at the bulk and all the current is now flowing from 
the incompressible bulk, slightly asymmetric with respect 
to the center. At the highest B field strengths shown in 
Fig|3] the system is out of the QHP and both the elec- 
tron and the current densities become symmetric, again. 
If compared to Ahlswede experiments the sequence of the 
potential types is I-I-III-III-II-I, which perfectly matches 
with the findings. It is reasonable to expect that, if one 



starts already with an asymmetric density profile, the 
asymmetry induced by the large current will be either 
enhanced or suppressed depending on the current direc- 
tion. A (+) bias will tilt the LLs resulting a high poten- 
tial on the RHS, whereas a (-) bias will do the opposite. 
Therefore, for the (+) bias the sequence is I-III-III-III- 
II-I, which essentially means that the asymmetric sam- 
ple enters to the QHP at a lower B field compared to a 
generic sample, due to already existing large IS at the 
RHS. The situation is rather different for the (-) bias, 
since the narrow IS is on the LHS and the high bias will 
enlarge this IS. Hence, there is a competition between 
the slope of induced Hall potential and the confinement 
potential to generate a wide IS. Thus, the asymmetric 
sample when (-) biased will enter to the QHP at a higher 
B field value compared to both (-I-) biased and generic 
samples. The asymmetric distribution of the ISs is obvi- 
ous for the generic sample, where we only show the ISs 
(dark-blue regions) wider than A. Our findings point out 
that the extend of the QHPs depend strongly on the sam- 
ple asymmetry and the current polarization. The experi- 
mental manifestation of the predicted rectification effect 
requires, first of all, high mobility (> 1.0 x 10^ cm^V/s) 
and narrow (2c? < 10 fim) asymmetrically depleted sam- 
ples. One possible option is to define the Hall bars similar 
to the ones investigated in Ref. |13j . at which, an asym- 
metry in Rl is observed. However, the effect is not pro- 
nounced to draw clear conclusions. The main drawback 
of the gate defined samples relies on the fact that, by 
using gates one cannot create very steep edge potential 
profiles, therefore, rectification is somewhat suppressed. 
Meanwhile, one can define Hall bars with steeper edge po- 
tentials by deep etching, of course, with different etching 
dept on both sides. However, it is known that etching can 
cause inhomogeneities at the density profile which may 
become important when considering narrow samples. A 
hybrid solution, i.e. one side etch, other side gate de- 
fined, seems to be the most reasonable solution. To ob- 
tain the extreme sharp edge on one side, it is desirable to 
perform the suggested experiments on cleaved edge over- 
grown (CEO) samples, where it has been shown that no 
ISs reside at the sharp edge [16] . The experiments need 
not to be done at very low temperatures (0.4 < T < 4.0 
K), whereas the imposed current should not exceed the 
breakdown current due to Joule heating [9 , which can 
easily be determined by the experiments. 

In conclusion, for the high mobility, narrow and asym- 
metric samples we predict that, the large current either 
enlarges or shrinks the QHPs depending on whether the 
asymmetry induced by the current and the asymmetry 
caused by the edge profile coincides or not. Based on 
our findings, we proposed three set of sample structures 
where the effect of the current induced asymmetry and 
thereby the rectification of the QHPs can be controUably 
measured. As a final remark, we note that at the edge 
IQHE regime, i.e. Type HI, a highly non-equilibrium 



situation is present, due to the competition between the 
enhancement of the ISs resulting from the laxge current 
and suppression due to steep potential profile, therefore 
we expect a hysteresis like behavior in this regime both 
depending on the sweep rate and direction of the B field 
and current amplitude. 
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